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Systems with long range interactions display some anomalies when its dynamics and thermo¬ 
dynamics are studied below certain conditions. Among these anomalies are the quasi-stationary 
states, which are exacerbated because of special initial conditions that are used here. We present 
in this letter a new Hamiltonian whose potential is inspired in the two-dipole interaction. An an¬ 
alytical solution is obtained for the equilibrium in the canonical ensemble that is coincident with 
the one obtained from computational simulations. However, results from this model presents a kind 
of nonequivalence of ensembles in long-living states before arriving to equilibrium. Thus, a com¬ 
plete characterization is made for the nonequilibrium through molecular dynamics. In which, novel 
quasi-stationary states are observed due to the long range interactions. 
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From years ago, the Ising model has been considered 
the most relevant tool to study magnetic properties and 
the statistical behavior of many-body systems. Extra 
efforts have been made to propose several variations of 
the Ising model, specially for theoretical and numerical 
modeling of systems with long range interactions. 

Typically, it has been accepted that at low energy in 
this kind of systems arises a phase identified by the pres¬ 
ence of a single cluster of particles floating in a diluted 
homogeneous background. At high energy a homoge¬ 
neous phase is recovered; the cluster disappears and the 
particles move almost freely. In the transition region, 
the system is characterized by the microcanonical en¬ 
semble with a negative specific heat: the resulting insta¬ 
bility is extremely relevant)] because of its strong impli¬ 
cations on experimental and theoretical features. This 
corresponds to an apparent thermodynamical inconsis¬ 
tency, which has been solved by Hertel and Thirring[T]. 
They proposed that the canonical and microcanonical 
ensembles are not equivalent close to the transition re¬ 
gion. Until now, this proposal has been successfully con¬ 
firmed by numerical simulations of systems with long 
range interactions[5HB]. 

Additionally, specialists have been developing various 
methods |16IH9) to manipulate in computational simu¬ 
lations and theoretical descriptions adequately the inter¬ 
actions in order to obtain a proper characterization of 
the behavior of these systems because of the wide range 
of applications that we find in nature that going from 
microscopic to macroscopic scale. For instance, the pos¬ 
sibility of controlling the matter at the molecular level to 
develop nano-machinery is a dream from decades ago. In 
the other limiting case, the possible understanding of the 
behavior in the astronomical scales is other motivation to 
catch deeply this kind of systems. Certainly the systems 
in nature are idealized without interactions, but a deeper 
understanding needs to include interactions, which can 
be extracted from this class of pivotal descriptions. 

Several mean held models have been recently intro¬ 
duced to study these anomalies related to long-range in¬ 
teractions. In part, the behavior of the kinetic energy 


and other thermodynamic observables are used to char¬ 
acterize the stationary states, where long-living states 
before arriving to equilibrium are observed |10H12j . One 
of these intriguing models is the Hamiltonian mean held 
model, whose properties become paradigmatic and per¬ 
tinent to characterize general systems with long range 
interactions. However, this behavior seems to be persis¬ 
tent and it has been observed in several self-gravitating 
models and other systems with long range interactions. 

The main goal in this work is to study anomalous be¬ 
havior of systems with long range interactions from a new 
Hamiltonian mean field model inspired in the two-dipole 
interactions. We introduce this model increasing the 
stock of mean held models. We expect to bring other non¬ 
standard properties and to correctly characterize them. 
The dependence on the orientations of the dipoles is char¬ 
acterized by the zeroth-order approximation of the dipo¬ 
lar interaction. In this work, we use the canonical en¬ 
semble to calculate the free energy, magnetization and 
internal energy for the equilibrium. Through numeri¬ 
cal simulation we characterize the mean kinetic energy, 
distributions, caloric curve, mean square of the displace¬ 
ment to obtain the diffusion law. We observe that due 
to out of equilibrium states, we can define two different 
quasi-stationary states (QSS). This anomaly is a novel 
observation in comparison with other models known in 
literature. 

Now, for theoretical and numerical modeling of sys¬ 
tems with long range interactions, we take a system of N 
identical coupled particles, with mass equal to 1, whose 
dynamics evolves in a periodic cell described by a ID 
Hamiltonian given by 

N 2 ^ 

H = YlY + ^Yl[cos{e,-9j)-3cos9iCos0j-Aij] ( 1 ) 

where Pi and 0i represent the momentum and the angle 
of orientation of the particle i, with i = 1, - ■ ■ N, being N 
the size of the system. The parameters e and Aij stand 
for the coupling and initial conditions. Meanwhile, the 
parameter is defined to consider properly the zero 
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of the energy as follows: 


Aij = cos ( 6 »oj - ) - 3 cos 6*oi cos 6>oj, (2) 

thus, the set of angles { 0 ofe} stands for the initial orienta¬ 
tions of the dipoles. The interaction coupling is rescaled 
by the number of particles to make the potential ther¬ 
modynamically stable[]. If e is positive, the system is 
ferromagnetic, but if e is negative the system is antiferro¬ 
magnetic. The equilibrium state can be exactly derived, 
however the QSS are not standard and cannot be exactly 
derived. Complementary to this, the implementation of 
numerical methods comes being an acceptable tool at the 
moment of studying and characterizing these anomalies. 
The spin vector related to each particle is given by 

r^i = (cos Si, sin 0i) (3) 

Therefore, we can introduce the total spin vector 

= {M^,My) = Mexp{i(l)) (4) 

i-l 

where {M^, My) and M are the components and the mod¬ 
ulus of the vector M, respectively, (fi stands for the phase 
of the order parameter. The equation of motion is 

Pi =-^ {2MxSm 9i +My cos 0i) ( 5 ) 

and the potential can be written as follows 

^=^( 2 M 2 -M 2 -A), ( 6 ) 

where A = In. the canonical ensemble, the 

partition function 

Z{I3, N)=J 6.% = Zk{)3, N)ZviP, N), ( 7 ) 

where Zk{I3,N) is the kinetic part of the integral and 
the Zy{I3,N) the interacting part. Therefore, 

Zk{I3,N) = J d%exp = (y) 

In addition, considering the Eqs.(|^ and Q we can write 

Zv{/3, N) = j 6.% exp - 2)^ (9) 

if { 6 * 0 ,i = 0} for all i, A — —2. Now, taking into ac¬ 
count, both the real and complex Hubbard-Stratonovich 
transformations [201 El], we obtain 

/ OO 

dxe-^'^^^\o{2/3ex)^ (10) 

-OO 

where Ikiy) is the modified Bessel function of kt/i-order. 
This last integral can be evaluated through the saddle 



FIG. 1: Equilibrium magnetization as a function of internal 
energy U. Triangles are data obtained from the simulation. 
The solid line is the analytical solution from canonical ensem¬ 
ble: U = ^ + 


point method in the thermodynamic limit, iV —>■ oo. The 
free energy per particle Lp is given by 

(/7(/3, N) = l^ In ;:^-|-e/3-|- inf[-/3ea;^-|-lnIo(2/3ea;)] (11) 

z 27r x>o 

The solution of the extremal is obtained by the transcen¬ 
dental equation. 


Ii(2/3ea:) 
Io(2^ea:)' 


( 12 ) 


The procedure is particular to this problem, but in part 
analogous to the Hamiltonian mean field model. The 
critical temperature is Tc = 1, that is the double com¬ 
pare with the Tc obtained for the Hamiltonian mean field 
model. The clustered phase is found for T < Tc and the 
homogeneous phase occurs for T > Tc- If e < 0, the 
equation has a trivial solution, cc = 0. By contrary, if 
e > 0, the equation has a set of values for x and /3, which 
defines the solution of the problem. Finally, it is ob¬ 
tained the internal energy per particle, as a function of 
the temperature and magnetization, as 


dp 


1 

w 


+ l-M^ 


(13) 


where M, is the solution of the extremal problem. 

Numerical simulations are carried out, by microcanon- 
ical molecular dynamics, to check the validity of the an¬ 
alytical results, by the canonical ensemble. Solving nu¬ 
merically the equations of motions considering the spe¬ 
cial initial condition called water bag initial conditions 
(WBIC), we show that the equilibrium is well described 
by the canonical ensemble. The magnetization M is the 
clustering degree of the particles. It is possible to define 
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FIG. 2: In the figure it is depicted the evolution of the kinetic 
energy, it is shown two QSS previous to the equilibrium. The 
first appears for 10 < t < 10® mean kinetic energy 2 < K > 
/Nqssi = 0.775 (dashed line). The second occurs for 10"^ < 
t < 4 ■ 10® with 2 < K > /Nqss2 = 0.761 (dashed-dotted 
line) and finally T^q = 0.948 (dotted line). The data are the 
mean value of 100 samples with N = 8192 and U = 1.38. 



FIG. 3: It is depicted the 2 < K > /Nqss2 for QSS 2 (cir¬ 
cles) and its correspondingly equilibrium values (triangles) as 
a function of the internal energy U. In consequence, equilib¬ 
rium values match the Te, (solid line), which coincide with 
the analytical solution in the canonical ensemble. Simulation 
takes N = 4096. 


the critical temperature Tc where the continuous trend 
of M vanishes. 


The magnetization M is obtained from the Eq.(13l as 
a function of the internal energy U, which corresponds 
to the solution we derive from the canonical ensemble. 
Correspondingly, we obtain solution by simulations in the 
microcanonical ensemble for several energies. In Figj^we 
depict the equilibrium magnetization M as a function of 
the internal energy U. Numerical data are represented by 
triangles. The analytical solution is depicted by solid line 
from canonical ensemble given by Eq. The critical 
point is located at Uc = 3/2 just twice the value obtained 
for the Hamiltonian mean field model. 


Other challenge to characterize the current model is to 
evaluate the kinetic energy. The behavior of this thermo¬ 
dynamic quantity is relevant to define properties of the 
system. In the dynamics of the model is it is important 
to observe where the kinetic energy is constant. As far as 
we know, this is relevant to search the behavior of other 
thermodynamic observables. In Figj^ we represent the 
dynamics of the system by the mean value of the kinetic 
energy in 100 samples with the size N = 8.192 particles. 
As said before, regions where the kinetic energy is stable 
are interesting for the present study. First (10 < t < 10®) 
and second (lO"* < t < 4 x 10®) regions define the called 
QSS. In this model we clearly identify two time inter¬ 
vals where the kinetic energy is constant, previous the 
last one where the equilibrium is reached. The value of 
first stable kinetic energy is upper than the second value. 
We use U = 1.38 because the system has the maximum 
anomaly when the stability of the kinetic energy is ob¬ 
served. For this value of energy U, the temperature that 


is related with the kinetic energy in the equilibrium, it is 
obtained analytically Teq = 0.948 and coincides with the 
numerical data with 0.3% of exactness. 

After observing the difference among the value of the 
kinetic energy in equilibrium with the corresponding val¬ 
ues out of equilibrium, these are the two QSS, we choose 
the deepest numerical result to compare with the equi¬ 
librium that we obtain from two ways, the theoretical 
canonical ensemble and the numerical microcanonical 
simulation. Therefore, in the Figj^ we superpose three 
curves that represent twice kinetic energy per particle 
2 < K > /Nqss 2 as a function of internal energy U. 
Triangles stand for equilibrium data (from the numeri¬ 
cal simulations) that coincide exactly with the analytical 
solution. But, circles correspond to the QSS in the sec¬ 
ond region. Here, we notice that there is a region where 
one set of values is not coincident with the other. In 
particular, there is a point, specifically U = 1.38 where 
the difference is the greatest. This point has been used 
to characterize the QSS. In the interval from 1.1 to 1.38 
we observe a disagreement between molecular dynamics 
simulation and canonical ensemble treatment in the limit 
o the mean field approximation. 

In order to fully characterize the dynamics of the sys¬ 
tem, we are interested in observing the behavior in phase 
space. Thus, for instance, we can take some pictures at 
several states of the evolution. Thus, in Figj^it is shown 
some snapshots at different times to characterize the dis¬ 
tribution in phase space; namely, t = 20 close to WBIC, 
at t = 200 the first region of QSS, at t = 600.000 the 
second region of QSS and at t = 5.000.000 the equilib¬ 
rium. In Figj^ (a), we see the distribution slightly de¬ 
formed from the WBIC with defined regions occupying 
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FIG. 4: Snapshots of the distribution function for the dipolar 
Hamiltonian mean field model into phase space, (a) For t = 
20 weakly deviated from WHIG, (b) For t = 200, the system is 
in the first region of QSS. (c) For t = 600.000, the system is in 
the second region of QSS. Finaly (d) For t = 5.000.000, the 
system reaches the equilibrium where ellipses characterized 
trajectories into phase space. 

the phase space. In FigQft), the distribution is extended 
by all phase space in a irregular form. In FigQc), the 
distribution takes paths enough defined showing certain 
regularities that we see in the representation into phase 
space. In addition, we see in Fig[^d) trajectories that are 
defined by elliptical shapes represented by several colors. 

In addition, in Figj^ it is represented a typical equi¬ 
librium distribution in phase space, obtained from the 
average of among various snapshots in similar conditions 
that we show in FigQd). 


FIG. 6: The dynamics is illustrated by the evolution of the 
variance. Anomalous diffusion is obtained in two regions that 
coincide with the two QSS previously discussed. The equilib¬ 
rium is obtained when 7 = 1. 

Finally, we represent the dynamics of particles of the 
system, by the variance of the angular displacement de¬ 
fined as crg{t) = — 9i{0)y/N. Equations of mo¬ 

tion are obtained from the Eq. ([^ that are numerically in¬ 
tegrated. We choose a combination of parameters where 
it is possible to observe dynamics at initial conditions 
out of equilibrium and that progressively acquires states 
with anomalous diffusion until the equilibrium is reached. 
Specifically, two superdiffusive states are observe in Fig|^ 
where ag oc t'*' and 7 > 1 . 
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FIG. 5: Equilibrium distribution is depicted in phase space. 
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